Photoemission Quasi-Particle Spectra of Sr 2 Ru04 
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Multi-band quasi-particle calculations based on perturbation theory and dynamical mean field 
methods show that the creation of a photoemission hole state in Sr2Ru04 is associated with a highly 
anisotropic self-energy. Since the narrow Ru-derived d XZtVZ bands are more strongly distorted by 
Coulomb correlations than the wide d xy band, charge is partially transferred from d xz<yz to d xy , 
thereby shifting the d xy van Hove singularity close to the Fermi level. 

PACS numbers: 79.60.Bm, 73.20.Dx, 74.70. Ad, 74.80.Dm 



Angle-resolved photoemission spectroscopy is one of 
the key techniques providing detailed information on the 
Fermi surface topology of high temperature superconduc- 
tors. In the layered copper oxide compounds (cuprates) 
high-resolution photoemission spectra can be obtained 
below and above the superconducting transition temper- 
ature and for different hole doping regimes. Although de 
Haas-van Alphen (dHvA) experiments in principle yield 
more reliable bulk Fermi surface data, they are less use- 
ful for the investigation of cuprates since they require 
extremely pure samples. Thus, it has so far not been 
possible for any of the high T c cuprates to obtain con- 
sistent Fermi surface data from both photoemission and 
dHvA measurements. 

The detection of superconductivity in Sr 2 Ru0 4 § is 
of great importance since this system is the only lay- 
ered perovskite compound known so far that is super- 
conducting in the absence of copper and without requir- 
ing doping. Thus, a critical comparison of photoemission 
Fermi surface data with those derived from dHvA mea- 
surements is feasable. Surprisingly, independent studies 
of the dHvA effect ||] and angle-resolved photoemission 
[^-|| yield highly contradictory Fermi surface topologies. 
This discrepancy raises serious questions concerning the 
interpretation of photoemission data also in cuprate su- 
perconductors. 

Because of the layered structure of Sr 2 Ru04, the elec- 
tronic bands close to the Fermi level may be qualitatively 
understood in terms of a simple tight-binding picture. 
These bands are derived mainly from Ru ti g states. The 
wide xy band exhibits two-dimensional character, while 
the narrow xz and yz bands are nearly one-dimensional. 
All three bands are roughly 2/3 occupied, giving about 
4 Ru d electrons per formula unit. Density functional 
calculations based on the local density approximation 
(LDA) [UQ| place the (ir, 0), (0, n) saddle point van Hove 
singularity of the xy band about 60 meV above the Fermi 
energy. Taking into account gradient corrections slightly 
lowers this singularity to about 50 meV above Ep 0. 
Fig. 1 provides a qualitative picture of the t2 g bands and 
of the Fermi surface exhibiting one hole sheet (a) and 
two electron sheets (/3, 7). Whereas the dHvA data ||] 
are consistent with these results, photoemission spectra 



reveal a fundamentally different topology the xy 

van Hove singularity near M appears below the Fermi 
level, converting the 7 sheet from electron-like to hole- 
like. Nevertheless, both experiments are reported to be 
in accord with Luttinger's theorem. 
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FIG. 1. (a) Dispersion of t zy bands of S^RuCU in simpli- 
fied two-dimensional Brillouin Zone (Ef = 0). (b) Solid lines: 
Fermi surface consistent with LDA band structure and dHvA 
measurements, with hole sheet a and electron sheets /3, 7 (af- 
ter accounting for hybridization). Dashed line: approximate 
xy Fermi surface derived from photoemission, indicating that 
7 is hole-like. 

Various effects have been proposed to explain this 
discrepancy between dHvA and photoemission data, 
(a) Since photoemission is surface sensitive, a surface- 
induced reconstruction associated with the breaking of 
bonds could lead to a modification of Ru-0 hopping and 
a shift of tig bands. Estimates based on bulk calculations 
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indicate that this mechanism could push the xy van Hove 
singularity below Ep H . Actual surface electronic struc- 
ture calculations for Sr 2 Ru04 ||, however, show that 
this singularity remains above Ep even in the presence 
of surface relaxation, (b) Slightly different doping levels, 
temperatures, and magnetic fields used in both experi- 
ments could result in different Fermi surfaces, but these 
effects are believed to be too small to explain the main 
discrepancy (c) That the interpretation of the dHvA 
data is incorrect seems unlikely in view of the large body 
of experience available for this method, (d) The interpre- 
tation of the photocmission spectra, on the other hand, 
is non-trivial because of the quasi-particle character of 
the hole state created in the experiment. This aspect of 
the data has been ignored so far and is the focus of the 
present work. 

According to the reduced dimensionality of Sr2Ru04, 
creation of a photohole should be associated with highly 
anisotropic screening processes which reflect the nature 
of the different electronic states involved. As shown in 
Fig. f, the relevant bands near Ep comprise a roughly 
3.5 eV wide band formed by in-plane hopping between 
Ru d xy and O 2p orbitals, and 1.4 eV narrow d xz , d yz 
bands. Assuming an on-site Ru dd Coulomb interac- 
tion U ~ 1.5 eV, we have the intriguing situation: 
W XZ} y Z < U < W xy , where Wi is the width of the i th 
t2 g band. A value U «s 1.5 eV was in fact deduced 
from the observation of a valence band satellite in res- 
onant photoemission from Sr2Ru04 |1C[ |, According to 
this picture, intra-atomic correlations have a much larger 
effect on the xz, yz bands than on the xy band, giving 
rise to a strongly anisotropic self-energy. Because of the 
~ 2/3 filling of the xz,yz bands, their narrowing, com- 
bined with Luttinger's theorem, leads to a charge flow 
from the xz,yz bands to the xy band. As we discuss 
below, for reasonable values of U this charge transfer is 
large enough to push the xy van Hove singularity close 
to or even below the Fermi level. 

Since we are concerned with the qualitative influence of 
multi-band correlations on quasi-particle spectra, we con- 
sider for simplicity next-nearest-neighbor tight-binding 
bands of the form e{k) = — EQ — 2t x cos ak x — 2t y cos ak y + 
At' cos ak x cos ak y , where (eqi t x , t y , t') = (0.50, 0.44, 
0.44, -0.14), (0.24, 0.31, 0.045, 0.01), (0.24, 0.045, 0.31, 
0.01) eV for xy,xz,yz, respectively (see Fig. 1). These 
parameters ensure that the xy band has edges at —2.8 
and 0.7 eV, with a van Hove singularity at 0.05 eV, and 
the xz,yz bands have edges at —0.9 and 0.5 eV, with 
van Hove singularities at —0.80 and 0.26 eV, in agree- 
ment with the LDA band structure 0. 

Next we specify the on-site Coulomb and exchange in- 
tegrals which we use in the self-energy calculations dis- 
cussed below. In the present case involving only t 2g 
states, there are three independent elements (i ^ j) [p"i"|| : 
U = U' = and J = (ij\\ji) = (H\\jj) = 

(U—U')/2, where i = 1 . . . 3 denotes xy, xz, yz. Thus, the 
Hartree-Fock energies are S^ F = n\XJ + 2n 2 (2U' — J) 
and T,fJ = m(2J7' - J) + n 2 (U + 2U' - J). As the 



band occupations rii are rather similar, it is convenient 
to define the average occupation n, so that n\ = h — 25, 
n 2 3 = n + S, and Sf F = 5n(U - 2J) + 2S(U - 5J), 
E£f = 5n (jj „ 2 j) _ S(U - 5J). 

It is instructive to consider the second-order contribu- 
tion to the local self-energy since the key point, namely, 
the large difference between the quasi-particle shifts of 
the xy and xz, yz bands, can already be illustrated in this 
approximation. Because the t 2g bands do not hybridize, 
the self-energy has no off-diagonal elements. The imag- 
inary parts of the diagonal second-order Coulomb and 
exchange terms are given by 

ImE.H = nJ^RjM (v\\ kl ) [2<W||*i> - (W||j»>] (1) 
jkl 

where 

(roo pO i>0 />0 poo roc 

/ / / + / / / ) duj 1 duj 2 duj3 
JO J-ooJ-oo J -co JO Jo ' 

x pj^^pk^^pi^z) 8(uj + ujx -u>2~ w 3 ). (2) 

Here, Pj(w) denotes the single-particle density of t 2g 
states. Exploiting the symmetry properties of the 
Coulomb matrix elements, Eq. ([!]) reduces to 

ImSiH = U 2 R nl {uj) + 2J 2 R l22 {uo) 

+ 4(C/' 2 + J 2 - U' J) ifciaM (3) 
Im Ea, 3 (w) = {U 2 + 2U' 2 + 3 J 2 - 2U'J) R 222 (oj) 
+ J 2 i?2iiH 

+ 2(U' 2 + J 2 - U 1 J) R ll2 {uo) . (4) 

The above expressions demonstrate that even for J = 
the self-energy of a given band depends on scattering 
proccesses involving all three t 2g bands. Nevertheless, 
S xa is dominated by interactions within the wide xy band 
or between xy and xz, yz. On the other hand, Y* xz , yz pri- 
marily depends on interactions within the narrow xz, yz 
bands or between xz, yz and xy. These differences are 
a consequence of the layered structure of Si'2Ru04 and 
give rise to anisotropic relaxation shifts. 

For a more accurate description of charge transfer 
among quasi-particle bands, we include self-consistency 
in the spirit of dynamical mean-field theory |l2[ . In this 
scheme, Si is a functional of the effective bath Green's 
function Q~ l — G~ x -I- Si, where the local d is given by 

Gi(w) = duj' ^9-1 7 ■ (5) 

A typical frequency variation of Sj is shown in Fig. 2. 
Near Ep, the imaginary parts vary quadratically with 
frequency and the real parts satisfy S a;Zj j )2 3> S xy , i.e., 
the energy shift of the narrow xz, yz bands is much larger 
than for the wide xy band. Moreover, the difference 
^xz,yz — S^j, at Ep is much larger than the difference 
between the Hartree-Fock energies ^ x ^ yz — ^ xy F ■ 
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Qualitatively similar results are derived from more re- 
fined treatments of on-site Coulomb correlations using 
multi-band self-consistent Quantum Monte Carlo (QMC) 
methods [|l2],|l3| . The temperature of the simulation was 
15 meV with 128 imaginary time slices and ~ 300 000 
Monte Carlo sweeps. Fig. 3 shows the quasi-particle den- 
sity of states Ni(ui) = — ^ Im Gj(a>), obtained via maxi- 
mum entropy reconstruction fL4| , together with the bare 
density of states Pi(uj). The van Hove singularities near 
the edges of the xz,yz bands arc shifted towards Ep, 
causing a sizeable band narrowing. Because of the ^2/3 
filling of these bands, this effect is not symmetric, giving 
a stronger relaxation shift of the occupied bands than for 
the unoccupied bands. There is also some band narrow- 
ing of the xy bands, but since U < W xy this effect is 
much smaller than for the xz, yz bands. 
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FIG. 2. Real and imaginary parts of self-consistent sec- 
ond-order self-energy for U — 1.2 eV, J = 0.2 eV. Solid 
curves: xy, dashed curves: xz. 

A crucial point is now that in order to satisfy the Lut- 
tinger theorem the more pronounced band narrowing of 
the xz, yz bands requires a transfer of spectral weight 
to the xy bands. Thus, the xy van Hove singularity is 
pushed towards the Fermi level. In the example shown 
in Fig. 3, it lies about 10 meV above Ep, compared to 
50 meV in the single-particle spectrum. We emphasize 
that this result is a genuine multi-band effect where the 
filling of a relatively wide quasi-particle band is modi- 
fied by correlations within other narrow bands of a dif- 
ferent symmetry. Since the values of U and J are not 
well known, and considering the approximate nature of 
our single-particle bands and self-energy calculations, it 
is not possible at present to predict the exact position of 
the xy singularity. It is conceivable, therefore, that this 
saddle point might lie even closer to or below Ep. 

As indicated in Fig. 1, the topology of the Fermi sur- 
face of Sr2Ru04 depends critically on the position of the 
xy van Hove singularity with respect to Ep. It is ev- 
ident therefore that the charge transfer from xz, yz to 
xy due to the creation of the photohole must be taken 
into account when using angle-resolved photoemission 



to determine the shape of the Fermi surface. To com- 
pare our results with photoemission spectra, we show in 
Fig. 4(a) the dispersion of the t 2g quasi-particle bands 
along TM and MX derived from the spectral function 
Ai(k, u) = — ~ Im [uj + p — £j(k) — £i(a>)] . The xy van 
Hove singularity at M lies 10 meV above Ep, so that 
considerable spectral weight appears below Ep in the im- 
mediate vicinity of M. To account for the finite energy 
resolution, and following the experimental procedure for 
determining the spectral weight near Ep we show in 
Fig. 4(b) the Fermi surface obtained from the partially in- 
tegrated spectral function A,(k) = f A A du> Ai(h, uj + iA) 
with A = 25 meV. Considering in addition the finite 
aperture of the detector (typically ±1°, corresponding 
to ±5% of fcii near M for 25 eV photon energy), it is 
unavoidable to pick up spectral weight from occupied re- 
gions near M, even when the detector is nominally set 
at M. Thus, the near-degeneracy of the xy singularity 
with Ep makes it extremely difficult using angle-resolved 
photoemission to determine the fc-point at which the xy 
band crosses the Fermi energy. Photoemission data taken 
with better energy and angle resolution might provide a 
more conclusive answer. 
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FIG. 3. (a) Quasi-particle density of states Ni(uj) derived 
from self-consistent QMC scheme for U — 1.2 eV, J = 0.2 eV. 
(b) Single-particle density of states pi(u)) derived from tight 
binding bands. Solid curves: xy, dashed curves: xz. 

Figs. 3 and 4 also show that due to the narrowing of the 
xz, yz bands, the weakly dispersive band is shifted from 
—0.8 eV to about —0.4 eV, in agreement with photoemis- 
sion data [||-||. For fcy between M and X, this band is 
observed to cross Ep at about (it, 0.67r), in good accord 
with our calculations. In addition, the calculations indi- 
cate the existence of a satellite below the xz, yz bands 
which might be related to the spectral feature observed 
near 2.5 eV binding energy using resonant photoemission 
fllPf . The precise location of this satellite is difficult to 
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determine because of the uncertainty of U and the ap- 
proximate nature of our self-energy calculations. 
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FIG. 4. (a) Quasi-particle bands along FM and MX de- 
rived from self-consistent second-order self-energy. Symbols: 
tight-binding bands, (b) Quasi-particle Fermi surface after 
accounting for energy broadening and resolution (see text). 

Because of the proximity of the quasi-particle xy van 
Hove critical point to the Fermi level, the imaginary part 
of the self-consistent self-energy exhibits a small linear 
contribution near Ep, indicating that the system may 
partially behave like a marginal Fermi liquid. In fact, in 
Eq. (3), it is only the first term ~ Rm((J) that gives rise 
to a linear term if the singularity coincides with Ep. As a 
result of multi-band effects, however, this contribution is 
rapidly dominated by stronger quadratic terms involving 
the narrow xz, yz bands. Thus, we find the marginality 
to be rather weak. 

We finally discuss the mass renormalization derived 
from our quasi-particle bands. For Coulomb and ex- 
change matrix elements in the range U = 1.2 — 1.5 eV, 
J = 0.2 — 0.4 eV we find ra* jm w 2.1 — 2.6, in agreement 
with photoemission estimates m*/m w 2.5 B. 
dHvA measurements § and specific heat data 
gest a factor of 3 — 4. 

In summary, multi-band quasi-particle calculations for 
Si'2Ru04 show that the simultaneous existence of nearly 
one- and two-dimensional t2 g bands near Ep leads to a 
highly anisotropic self-energy of the photoemission hole 




state. Because of Luttinger's theorem, this anisotropy 
gives rise to a charge flow from the narrow xz, yz bands 
to the wide xy band, thereby shifting the xy van Hove 
singularity very close to Ep. As a result, in the vicin- 
ity of M considerable spectral weight appears below Ep. 
These results might explain the controversial nature of 
recent photoemission data which have difficulty in deter- 
mining whether or not the xy band at M is occupied. 

The calculations were performed on the Cray T3e of 
the Forschungszentrum Jiilich with grants from the John 
von Neumann Institute for Computing. 
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